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O ■ Abstract 

In the model III 2HDM there are new CP violating phases which would affect the s — B® g 
mixing. In this paper, we calculate the new physics contributions to the neutral B meson mass 
splitting AMg (q=d, s) at the next-to-leading order (NLO) level. Using the high accuracy data 
and other relevant data, we draw the constraints on the parameter space of the model III 2HDM. 
Moreover, we calculate the new physics corrections to the ratio q/p. It is found that the phase of 
■ (q/p)n for Bd which is due to the new contributions is very small and consequently in agreement 

with the measurements of the time dependent CP asymmetry Sj/^k hi B — > J/ipKs- On the 
| contrary, the phase of (q/p) n for B s is large enough to give significant effects on CP violation in 

^ ■ the neutral B s system. 
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I. INTRODUCTION 



The neutral B-meson system has been of fundamental importance in testing the Stan- 
dard Model (SM) picture of flavor changing neutral currents (FCNC) and CP violation 
and in probing virtual effects from potential new physics beyond SM at low energies. The 
presentence of CP violation in the neutral B meson system has been established. The 
measured value |2j,|3j of the time dependent CP asymmetry Sj/^k in B — > J/ipKs 

S mK = sin(2/5( J/^))woHd-avc = 0.734 ± 0.054. (1) 

is in agreement with the prediction in the standard model (SM). 
In the B® (q=d, s) system, we expect model independently that 

r 12 < m X2 (2) 



of B® — Bg. Thus one has 



where Ti2 and M\ 2 are the off-diagonal terms of 2 x 2 decay width matrix and mass matrix 

(3) 



q M* 2 



where q and p are defined as 



P |M 12 |' 
AM Bq = 2\M 12 \, (4) 



\B l ,h) =p\B°)±q\B°) (5) 



with Bi h denoting the light and heavy meson mass eigenstates, and the normalization 
condition is 

\p\ 2 + \q\ 2 = l. (6) 
In the SM, one has, according to the box diagram calculation, 

1 = -Y^k ( 7 ) 

v v tb v t * q { ] 

and the deviation of \q/p\ from 1 is ~ 10~ 3 (10~ 5 ) for Bj (B s ) which is unobservably 
small. Eq. (JJj) for the B® system has been verified by measurements of the time dependent 
CP asymmetry Sj/^k in B — > J/ipK s , as mentioned above. Therefore, it would give a 
constraint on parameters of new theoretical models. 

As it is obvious from Eq. (@J), B° — B° mixing is responsible for the small mass 
differences between the mass eigenstates of neutral B mesons. The results of measurements 
of the mass splitting Am^ have been obtained with high accuracy. The current world 
averages of Am, (q = d, s) follow @,0 

AM Bd = 0.502 ±0.007 ps- 1 

AM Bs > 14.4 ps" 1 (8) 

which would impose a stringent constraint on new model building. 

In the SM, B° — B° mixing is dominated by the box diagrams with two internal t- 
quarks and W gauge bosons. In new physics models, the box diagrams with one or two 
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W gauge bosons replaced by the new charged scalars or vector bosons and/or top quarks 
replaced by new fermions can also contribute to B° — B° mixing. Using the precision data 
we can put stringent constrains on the parameter space of a new physics model. In the 
paper we study -B — -B mixing in the model III 2HDM. 

It is well-known that in the model III 2HDM the couplings involving Higgs bosons and 
fermions can have complex phases, the new CP violating phases would affect B° — B° 
mixing. B°-B° mixing has extensively been studied in the model III 2HDM during the 
past years. The charged-Higgs boson contributions to B°—B° mixing have been calculated 
at the leading order for a long time In the framework of SM, Ref. [9 J is the first one to 
present the NLO QCD corrections to B° — B° mixing while in conventional model I and 
model II 2HDM Ref. is. The possible constraints on the model III 2HDM from the 
measured parameter Xd = AMr/Tr were studied, for example, in Refs. HUH EDI at the 



LO level. At the NLO level, Ref. [1J| has calculated the charged Higgs boson contribution 
to mass splitting AM# d and drawn the constraints on the parameters of model III in terms 
of the high precision data, considering the uncertainty of the non-perturbative parameter 



fB d yBB d - The effects of the new CP violating phases on AM# d have been taken into 
account 13] . However, in Ref. [l3| as well as Ref. [lij all new parameters except for \tt,w> 
in the model III 2HDM are set to be zero. Furthermore, the ratio q/p which is important 
ingredient to study CP violation in the neutral B system has not been analyzed. In this 
paper, we will calculate the new physics contributions to AMb q taking into account the 
effects of the complex phases and keeping all relevant parameters of the model, which are 
only subjective to constraints from experiments, non zero at the NLO level in the model 
III 2HDM. By comparing the theoretical predictions with the high accuracy data, we draw 
the constraints on the free parameters of model III. Moreover, by using the constrained 
parameters we calculate the contributions of new physics to the ratio q/p in the neutral 
B system. 

The organization of this paper is as follows. In the next section, we describe the model 
III 2HDM briefly. In Section III, we give the effective Hamiltonian responsible for B° — B° 
mixing and calculate the mass splitting AMg ? and the ratio q/p at the NLO level in the 
model III. The section IV is devoted to numerical results. Finally, we conclude in section 
V. 



II. THE MODEL III TWO HIGGS DOUBLET MODEL 

As the simplest extension of the SM, the so-called two- Higgs- doublet models [u| 
may naturally have flavor changing neutral currents (FCNC's) mediated by the Higgs 
bosons at the tree-level, unless an ad hoc discrete symmetry is imposed. In this paper, 
we shall focus on the model III 2HDM Q UtI Il8^ . In the model III, there is no discrete 



symmetry and both the Higgs doublets can couple to the up-type and down-type quarks. 
In general one can have a Yukawa Lagrangian of the form 

C Y = VijQi.AU^n + n^Qu.Ih P.,.,* + C'jQ.j.HJ-jm + ^Qi, L H 2 D j>R + h.c. (9) 

where Hi (i = 1,2) are the two Higgs doublets, while rjff and (h3 = 1)2,3 are 
family index) are the nondiagonal matrices of the Yukawa couplings. We can choose to 
express Hi and H2 in a suitable basis such that only the rf^- D couplings generate the 
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fermion masses, i.e. 



V2 

Then the two Higgs doublets in the basis are of the form 



(# 2 > = (10) 



\ ,(,/2G* 
v + <t>l J V ^° 



^ 2 = 7f(0(l^ o) 



where iT are the physical charged-Higgs bosons and A is the physical CP-odd neutral 
Higgs boson, G° and G ± are the Goldstone bosons that would be eaten away in the 
Higgs mechanism to become the longitudinal components of the weak gauge bosons. The 
advantage of using the basis is that the first doublet Hi corresponds to the scalar doublet 
of the SM while the new Higgs fields arise from the second doublet H2. The 0° and 0° are 
not the neutral mass eigenstates but linear combinations of the CP-even neutral Higgs 
boson mass eigenstates, H° and h°: 

H° = 0° cos a + 0^ sin a (12) 
h° = — cj)\ sin a + 0° cos a 

where a is the mixing angle, such that for a = 0, (0°, 0°) coincide with the mass eigen- 
states. 

Note that Q^l, U^r and Dj R in Eq.® are weak eigenstates, which can be rotated into 
mass eigenstates. After the transformation the flavor changing (FC) part of the Yukawa 
Lagarangian becomes 3 



C y ,fc = Z',jQu. I hi).,; + ^Q l , L H 2 D 3tR + h.c. (13) 

where 

^D ={V U, D) - K ^D_ V U,D (M) 

in which V^' R are the rotation matrices acting on the up- and down-type quarks, with left 
or right chirality respectively, so that Vqkm = {Yl)^l ls the usual Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. Feynman rules of Yukawa couplings follows from Eq. (|T3*|) can 
be found in, e.g., Refs. jl^, EH- The FCNC couplings are given by the matrices ^ U,D and 
the charged FC couplings are given by 

^ charged > v ( 



charged ~ S ' V CKM 
■D 

charged 



-V C KM-i D (15) 

Because the definition of the £^' D couplings is arbitrary, we can take the rotated couplings 

as the original ones and shall write £ ,D in stead of £ U ' D hereafter. The Cheng-Sher ansatz 
for is 

by which the quark-mass hierarchy ensures that the FCNC within the first two generations 
are naturally suppressed by the small quark masses, while a larger freedom is allowed for 
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the FCNC involving the third generation. In the ansatz the residual degree of arbitrariness 
of the FC couplings is expressed through the \j parameters which are of order one and 
need to be determined by the available experiments, and they may also be complex. In 
the paper we choose ^ U ' D to be diagonal and set the masses of u and d quark to be zero 
for the sake of simplicity, so that Xu (i = c,s,t,b) are the new free parameters and will 
enter into the Wilson coefficients relevant to the process. 



III. B° - B° MIXING IN THE MODEL III 2HDM 



A. Effective Hamiltonian 



The most general effective Hamiltonian for AB = 2 processes beyond the SM can be 
written as 

5 3 

k a / =2 = E C *Q* + E w« (17) 



i=i i=i 



with 

Qi = <%rb a L <ihX 

Q2 = m.?A 

Qs = €A€ R bl 

q a = 

q 5 = rA€K (is) 

where q = s, d, corresponding to the operators of B s and B d system respectively, Pl,r = 
(1 =)= 7s)/2, and a, (3 are colour indexes. The tilde operators Qi {i = 1, 2, 3) are obtained 
from Qi correspondingly by exchanging L «-> R. 

At one loop level only the charged scalars are relevant for the box diagrams contributing 
to the B° — B° mixing amplitude. Using the Cheng-Sher ansatz for £ U,D which we assume 
in the paper, we rewrite the tree level couplings of = (H + , G + ) (where G ± is the 
would-be Goldstone boson) as follows 

C mt = H+u A V AI (Yf n P L + Yi n P R )dj + h.c. (19) 

where 

y aii = e m UA | \* AA for 1=1 

L ' V2s w M w \ 1 for 1=2 

V AII _ e m dr v / -A// for 1=1 

I D = — A 



R ' V2s w M w \ -1 for 1=2 

with A = t, c, u and I — b,s, d. 

In terms of the coefficients Y^ n and Y^ n we can express the contributions of to 
the Wilson coefficients Cj of the relevant operators responsible for B° — B° mixing as 
follows. 
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Diagrams with one W and one H give : 

yt^VAJl n / 2 ^,2 fl/f 2 ^2 \ 

C HMJ - 77~2 ^2~Z^ m u A Y L Y L U {m UA ,m UA ,M w ,M H+ ) 

1U7T Zijy ^ 

CbO*) = {VX ^ A f A E E 4 AJ ^4A)o « , < , M* w , < + ) (20) 

where svk = sin^ is the Weinberg angle), the definition of the four-point integral 
functions D and D 00 can be found in the Appendix. 
Diagrams with two Hf give 2 : 

fif..\_ 1 (^A/^Aj) 2 V y\AIl Y AJny\AIn Y AJl n / 2 2 ^2 ^2 s 

A l,n 

rjn\- 1 (^AJ^Aj) 2 V^V^ 2 ytAIZyAJnytATnyAJJ n /,„2 2 M 2 »^2 \ 
^2(^J--2 l^l^ m UA Y R Y L Y R Y L U {m UA ,m UA ,M H +,M H +) 

A l,n 

r n (,A- 1 ( V AI V Aj) 2 2 yfAUyAJnyt^riyAJi r, n ( m 2 m 2 ^2 ^2 s 

^2lMJ - -g — J^2 — 2^2^"Vt y £ y « y L Y R L, o{rn UA ,m UA ,M H +,M H +) 

A l,n 

rJn\- ( v ai v aj) 2 sr^sr^ 2 Y ^ AIl Y AJn Y^ AIn Y AJl r>J™ 2 ™ 2 M 2 M 2 \ 
2^2^ m «A y L Y L y r Y R 1J o{'m UA ,m UA ,M H +,M H +) 

A l,n 

CM = 2 {Va £ a / Y: f: Yl An Y A ^Y^Y AJl D 00 « , < , , ) (21) 
A l,n 

The above results of the Wilson Coefficients are in agreement with the ones in [22], the 
different sign and factor are due to the different definitions of the four-point integral 
functions. At one loop level there are no contributions to the Wilson coefficients of the 
operators Q 3 and Q 3 . 

In order to calculate the NLO QCD corrections to the Wilson coefficients Cj in the 
evolution from the scale of new physics fi t down to the low hadronic scale \x (~ one 
has to solve the corresponding renormalization group equations. The details can be found 
in Ref. [23j], with the result 

C r {n%**) = J2J2 + ^"O ^ °M ( 22 ) 

i s 

where we have set the new physics scale [it = mt and rj = a s (fj,)/a s (mt). The magic 
bf' s ^ and cf' S; 



numbers a«, b\ ' and c\ ' are given in Ref. |23[. 



1 The contribution of G ± to Ci(/i) is already taken into account in the Inami-Lim function So(xt) |2l| . 
Masses of the u and c quarks are neglected. 

2 In the sum over I and n in the expression for Ci (p) the contribution of G G^ is excluded since it has 
been taken into account in the function So(xt) [2l|. 
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B. The Master Formulae 



In terms of the matrix elements of the effective AB = 2 Hamiltonian we have the 
neutral B meson mass difference 



2|M- 



12 



AM Bq 

and the ratio of mixing parameters 3 

{q/p)B q = —exp \ —iArg[M 

with 



(?)■ 

12 . 



M- 



(?) 

12 



(B°\H eff (AB = 2)\B° q ). 



In the SM, the mass splitting AM» at the NLO level is available [2 



AM 



B a 



c 2 m 2 

^^ni Bq (B B jl q ) VB S (xm q V t l\ 2 



(23) 

(24) 
(25) 

(26) 



where Sq(x) is the Inami-Lim function |2l| . the details of B Bq and r\ B can be found in 
Ref. 0. 

In the framework of model III 2HDM, we need to know the matrix elements of the 
relevant operators Qi and Qi between neutral B mesons. One defines jii 



(B q \QM\B q ) 
(B q \Q 3 (ji)\B q ) 

(Bq\Qi(jl)\B q ) 

(B q \QM\B q ) 



1 

g ■ • -J->q J O q 

5 



mBjlB?^) 



m B a 



24 \m b (n) + m q (ji) 



24 \m b ((j) + m q (fjL) 
m Bq 

4 \m b (p) + m q (n) 
1 ( m Ba 



mBjlBf^) 



(«), 



mBjlBfin) 



(<?), 



mBj B B?(») 



(<z), 



(27) 



12 \m b ([i) + rrig(n) 

where f Bq is the decay constant and Bi is the so-called bag factor. The matrix elements 

of Qi-i are the same as that of Qi-3- The lattice calculations of B^ 9 \fj.) have been done 
in Ref. 12511 and the results are 



5{ d) K) = 0.87(4) +5 B [ 



Bf\m b ) = 0.82(3) (4) 
B^lmb) = 1.02(6) (9) ^ 



5l d) (m fe ) = 1.16(3)_ 7 

h2: 
-7 



+5 



5f (m b ) = 1.91(4)+ 22 



b { : 

b[ s 



(m 6 ) = 0.86(2)1* 
(m b ) = 0.83(2) (4) 
(m 6 ) = 1.03(4)(9) 
(m 6 ) = 1.17(2)1 5 7 
(m b ) = 1.94(3)lf 



(2f 



For the sake of simplicity we shall suppress the subscript B q of q/p hereafter. 
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IV. NUMERICAL ANALYSIS 



In numerical calculations, the following values of parameters are assumed: 

m Bd = 5.279 GeV, m Bs = 5.370 GeV, m H ± = 200 GeV (29) 

The values of f Bq and B Bq in the lattice QCD calculation and in QCD sum rules have 
been estimated in so many works, see, e.g. [1^, 22, HI, H^. Here, we quote the results 
listed in @. 

f Bd = i91±23l? 9 MeV 
f B . = 220 ±25 MeV 

f Bd \f^B~ d = 221±28t° 22 MeV 



fs s y B Bs = 255±31MeK (30) 

In [3| we calculated the constraint on |A&&| and \\ u \ due to the experimental upper 
bound of Br(B s — > H + which is p 

B r (B s -> /i+/i~) < 2.0 x 10~ 6 (CL = 90%) (31) 

Recently, the Ref. gave the new bound of Br(B s — > /i + /i~) which reads 

B r (B s -> fj+fj,-) < 5.8 x 10- 7 (CL = 90%) (32) 

With the updating upper bound of Br(B s — > /i + /i~), we recalculate the constraint on 
\Xbb\ and | Ah | and the result is shown in Fig.l. The cyan region is allowed by the new 
experimental bound of Br(B s — > Comparing with the Fig. 3 in Ref. [19], we find 

that the new bound gives a more stringent constraint on the values of |A&&| and \X tt \, as 
expected. 

It is shown in Ref. 0] that the strictest constraints on Abb and Xu come from B — > X s ^ 
and the neutron electric dipole moment (NEDM). Considering the theoretical uncertain- 
ties, we take 2.0 x 10~ 4 < By(B — > X s ^) < 4.5 x 10~ 4 , as generally analyzed in literatures. 
The NEDM can be expressed as 



<*(/*) J \0.lJ \M k± 



di = 10- 25 e-cm Im(A tt Abb) ( J ( 777 ) # ( ThH (33) 



with 



H(y) = * , V - fy-3-^y] (34) 

Here the parameter £ g has been estimated and has different values obtained by different 
methods: 0.03, 0.07 and 1, which is due to the hadronic uncertainty [l3|. Using the result 
in Ref. [3], considering the uncertainty of £ 9 and finding the maximally possible |AbbA tt | 
which satisfies the experimental constraints, and requiring that the charged Higgs boson 
is not too heavy (say, m#± < 250GeV) and the phase of XbbXu is not too limited, it follows 
that |A b bA tt | ~ 3. 
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From the well measured physical observable AM^ and using the values of |A&&| and 
| Aft | analyzed above, one can get the constraint on |A CC | in the model III. In Fig. 2, we show 
the correlation between |A CC | and 6^ {6u is the phase of Xu) for Qbb + Qtt — vr/2 4 due to the 
constraint from AM^ d within the la deviation (i.e., AMg^ = 0.502 ± 0.007ps _1 ). One 
can see from the figure that the constraint on |A CC | is very stringent and the dependence 
of |A CC | on 6bb is weak. We have also carried out the calculations for different values of 
Obb + Ou in the allowed range and the results are similar, i.e., the dependence of |A CC | on 
both 9bb and 6 U is weak. 

With the experimental lower bound of AMg s , i.e., AMg s > 14.4p.s~ 1 , we can draw the 
lower bound on |A SS |. Fig. 3 is the contour plot of |A SS | versus 9 SS for fixed \Xtt\, \Xbb\ and 
\X CC \ which satisfied the constraints from Bd — Bd mixing , T(b — > 37), T(b — > cru T ), po, Rb 
and the electric dipole moment of neutron. The yellow area is excluded by the measured 
experimental bound. Above the yellow region, all the values of |A SS | are allowed. We 
hope that the future B factory can make precise measurements of AMb s , consequently, 
the more stringent constraints on the free parameters in the model III can be obtained. 

Using the values of parameters obtained by analyzing above, we calculate the new 
contribution to the ratio q/p. The Figs. 4 and 5 are devoted to Arg[(q/p) n ] which denotes 
the new contribution to the phase of q/p versus 8 CC for B s and Bd systems respectively. 
As for the B s system, in Fig. 4, we show the Arg[(q/p) n ] dependence on 6 CC for fixed 
|A CC | = 100 and for |A SS | = 80 (solid curve), 100 (short-dashed curve) and 120 (dotted 
curve) respectively. From the figure we find Arg[(q/p) n ] increases with |A SS | increasing. 
In particular, Arg[(q/p) n ] is large enough to give significant effects on CP violation in the 
neutral B s system. In Fig. 5, for the Bd system, the solid curve corresponds \X CC \ = 98 and 
the dotted one corresponds lA^I = 101. We find the phase of (q/p) n is very small, which 
verifies the expectation in Ref . |l9j , and consequently in agreement with the measurements 
of the time dependent CP asymmetry Sj/^k in B — > J/ipK S - In contrast with the case 
of B s , Arg[(q/p) n ] decreases while |A CC | increasing. 

V. CONCLUSIONS 

In summary, we have calculated the new physics contributions to the neutral B meson 
mass splitting AM# ? (q=d, s) at the NLO level in the model III 2HDM, while taking 
into account the new CP violating complex phases involved Higgs bosons and fermions 
coupling. By comparing the theoretical predictions with the high accuracy data as well 
as other relevant data, we have drawn the constraints on the free parameters of model III. 
Moreover, by using the constrained parameters we calculated the corrections to the ratio 
q/p. It is found that the phase of (q/p) n for Bd which is due to the new contributions is 
very small and consequently in agreement with the measurements of the time dependent 
CP asymmetry Sj/^k in B —> J/tpK s . On the contrary, the phase of (q/p) n for B s is 
large enough to give significant effects on CP violation in the neutral B s system. 



4 The terms relevant to 8 CC in the formula of AMs d have been neglected because they are very small 
compared with other terms. 
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Appendix 



The Passarino-Veltman one-loop four-point functions with zero external momenta are 
defined as 



/- 

J % 



d A q 



= D (m , m 1} m 2 , m 3 ) (35) 



in 2 (q 2 — ml)(q 2 — m\)(q 2 — m 2 )(q 2 — ml) 

d^q 'llt'li- n , , . 2 ,,,2 , , 2 ...2 , 



in 2 (q 2 — ml)(q 2 — m\)(q 2 — m 2 )(q 2 — ml 



where 



1 1 



A)(m > m i> m 3> m 3) = fdo(x,y,z) 

O TTlo 



l l 

12 mf- 



g lxl/ D 00 (m ,m 1 ,m 2 ,m 3 ) (36) 

(37) 
(38) 



X = 


4v = 


2 

raj 
m\ 


z = 


coto toto 


X = 


m 2 

— 2>S/ = 
™3 


2 

raj 

2 ' 

m 2 3 


z = 





with 



f dQ (x,y,z) = 6 
+ 



a; In a; 



+ 



ylny 



(1 - x)(x - y)(x - z) (1 - y)(y - x)(y - z) 
zlnz 



(1 -z)(z-a;)(z-s/) 



(39) 



fdao{x,y,z) = -3 



x 2 lnx 



t/ 2 In j/ 



+ 



(1 — x)(x — y)(x — z) 
z 2 lnz 



(i - y)(y - x)(y - z) 



(1 - z)(z - x)(z - y) 



(40) 
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FIG. 1: The constraint on |A&&| and |A^| due to the new experimental upper bound on Br(B s 
wh ich reads Br{B s -> ^+ fj,~) < 5.8 x KT 7 (CL = 90%). 
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FIG. 2: The constraint on |A CC | due to the measured AM^ within la deviation. 
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FIG. 3: The constraint on |A SS | due to the lower experimental bound on AMg s . 
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FIG. 4: Arg[(q/p) n ] for B s system versus the CP violating phase 9 CC , for |A CC | = 100. The solid 
curve stands for \X SS \ = 80, the short-dashed curve stands for |A SS | = 100 and the dotted curve 
stands for |A SS | = 120 respectively. 
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FIG. 5: Arg[(q/p) n ] for system versus CP violating phase 6 CC for |A CC | = 98 (the solid curve) 
and for |A CC | = 101 (the dotted curve) respectively. 
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